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We study spin dynamics in a two-dimensional electron gas with a pure gauge non-Abelian spin- 
orbit field, for which systems with balanced Rashba and Dresselhaus spin-orbit couplings, and the 
(llO)-axis grown GaAs quantum wells are typical examples. We demonstrate the duality of the 
spin evolution and the electron density dynamics in a system without spin-orbit coupling, which 
considerably simplifies and deepens the analysis of spin-dependent processes. This duality opens 
a venue for the understanding of this class of systems, highly interesting for their applications in 
spintronics, through known properties of the systems without spin-orbit coupling. 

PACS numbers: 72.25.-b 



The understanding of spin dynamics in a two- 
dimensional (2D) electron gas with spin-orbit (SO) in- 
teraction is highly important both for theoretical and 
applied spintronics, [J i| including the design of devices 
with controlled spin transport. In many physically in- 
teresting situations the SO coupling can be elegantly de- 
scribed as an effective non-Abelian vector potential. 0- 
17| There exists a class of systems, where the SO coupling 
corresponds to a pure gauge non-Abelian field. Therefore 
it can be gauged away and the behavior of a physical 
system should map to that of a system without SO cou- 
pling. In practice, there are at least two systems of this 
sort widely investigated from the applied point of view. 
These are 2D electrons with balanced Rashba and Dres- 
selhaus couplings, and the electron gas in the (llO)-axis 
grown GaAs quantum wells. 

The coupled spin-charge dynamics is commonly de- 
scribed using the diffusion approximation, where the rate 
of the spin precession is much less than the momentum 
scattering rate. [IH Gil Currently, high- mobility 2D struc- 
tures, where the time scale of momentum relaxation is 
longer than the spin rotation time.pol [2lj became avail- 
able. Since here spins can make several turns between 
collisions with impurities, the conventional Dyakonov- 
Perel mechanism 18 is not applicable, [22j and another 
type of analysis is required. In the present paper we solve 
this problem for systems with a pure gauge SO coupling, 
including quantum effects due to the weak localization. 
We show for this general class of systems the existence 
of a duality of observables allowing the spin dynamics 
to be fully mapped to the density dynamics in systems 
without SO coupling. As a result, several regimes, in- 
cluding magnetic field dependence of the spin dynamics, 
can easily be explored using a single formula. 

We represent the Hamiltonian of a 2D electron gas with 
SO coupling as follows 3, 13] (the system of units with 



h = 1 is employed) 

H = / d 2 P ^+ (idi + Aif V + W [*+, *] , (1) 
2m J 

where ^(p) is a spinor field operator, and the func- 
tional W [^ + , \&] contains all spin-independent contribu- 
tions, including the external potential, electron-electron 
interactions, and, possibly, the electron-phonon coupling. 
Here m is the electron effective mass, and Ai are 2x2 
matrix- valued components of a non-Abelian SU (2) gauge 
field describing the SO coupling. In the broad class of 
systems of interest, Ai is a pure gauge, that is it can be 
removed by a local SU(2) transformation. The general 
form of a pure gauge vector potential, 

Ai = moii (h • a) , (2) 

corresponds to the following SO Hamiltonian 

H so = a (h • a) (k • v) , (3) 

where er is a vector of Pauli matrices, a is the SO 
coupling constant, h is a three-component unit vec- 
tor for the SO field direction, v is a 2D vector in the 
(x, y) plane, and ai = avi. The two practically impor- 
tant systems described by the Hamiltonian of Eq. J3|) 
are: (i) the balanced Rashba-Dresselhaus system [23J, |24| 
with h = (±l,±l,0)/%/2, v = (±l,±l)/ v / 2 i and (ii) 
the (llO)-axis grown GaAs quantum well. [25l428{ where 
h = (0, 0, 1) , v = (1, 0) with the well axes chosen with 
respect to the crystal axes as x || [110], y || [001], and z || 
[110]. Both systems are expected to demonstrate highly 
anisotropic spin relaxation times with the spin compo- 
nent along the h-axis having a very low relaxation rate, 
arising only due to a spin-dependent disorder. [2^ Spin 
currents in the thermodynamical equilibrium state. [301] 
being common for 2D electron systems with SO coupling, 
are absent [l6[ in the structures described by the Hamil- 
tonian in Eq. ([3|). 
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A local SU (2) transformation, which gauges away the 
above type of SO coupling is ^f(p) = U^^ r (p), where 



U_4 = exp [ima (h ■ a) (p ■ v)] 



(4) 



The transformation Q renders invariant all spin- 
independent quantities, such as the charge and current 
densities, while the spin density transforms covariantly: 



S = itr{«rU^(S • <r)U A }. 



(5) 



When the SO coupling is gauged away the dynamics 
of the transformed spin density S(p,t) reduces to the 
spin dynamics in the electron gas without SO interac- 
tion, which significantly simplifies the analysis. Then, 
the physical spin density S(p,t) is restored by 



(6) 



to obtain measurable results. Here we follow this guide- 
line and show that this approach allows to describe all 
regimes of spin dynamics on the same footing. 

We consider a 2D electron gas with a SO coupling of 
Eq. (UJ and, initially, a uniform spin density S produced, 
for example, by a static magnetic field B. At t = 
the magnetic field is released and the spin relaxes due to 
a disorder potential and other interactions. To describe 
this process we first eliminate the SO by the gauge trans- 
formation of Eq. (|4]). Using Eq. (|5|) we find that the ini- 
tial uniform physical spin density is mapped to the spin 
texture S (p, 0) = Sy (p, 0) + Sj_ (p, 0) , where the term 



Si, (p,0) = h(S-h) 



(7) 



being parallel to h, is untouched by the transformation 
and remains uniform, while the orthogonal to h part 
transforms into the helix structure [3lL I32I ] 

S ± (p, 0) = [S - h (S • h)] cos(Q hx -p)-(Sxh) sin(Q hx -p), 

(8) 

where Qh x = 2mav is the helix wave vector. 

Since in the transformed system there is no SO cou- 
pling, the uniform part of the initial spin distribution, 
Eq. ([7]). is constant in time. A nontrivial dynamics oc- 
curs in the orthogonal channel due to a diffusional decay 
of the initial helix spin texture described by Eq. J5]): 



S* (p, t) = J (p-p\ t)S? (p\ 0) d 2 p\ (9) 

where V^ 1 ^ 2 (p — p' , t) is the exact spin diffusion Green's 
function of a 2D electron gas, which takes into ac- 
count the disorder, electron-electron and electron-phonon 
interactions. To proceed further we note that in a 
nonmagnetic system without SO coupling the spin dif- 
fusion Green's function is diagonal in spin subspace 
T>^ 102 (p,t) = 5 01 p 2 V(p,t). Hence Eq. © simplifies as 



where T>(q, t) is a Fourier component of the spin diffusion 
Green's function 



V(q,t) = J d 2 pe~ i(q p) P(p,t), 



(11) 



and we have taken into account that only the Fourier 
components of T>(p, t) with the modulus of the wave vec- 
tor q = Qhx contribute to the dynamics of the helix in 
Eq. ©. Since the time-dependent factor in Eq. (fTU|) 
is scalar, the transformation back to the physical spin, 
Eq. (|6|), simply reduces to removing tildas and the coor- 
dinate dependence in Eq. (fT0|) . Thus, we get the following 
exact result for the observable spin evolution 



S± (t) = S ± (0) 



2^ 



V(Q hK ,u>)e 



-icut 



(12) 



In Eq. (fT2j) we represented V(q, t) via the Fourier inte- 
gral because in the w-domain there is a simple expression 
of the spin diffusion Green's function T>(q,uj) in terms 
of the spin-spin correlator (the spin response function) 

r c] 

Xpp{q,w), where (3 = (x, y, z) is the Cartesian index cor- 
responding to the spin component 



V(Q hx ,uj) = — 

1UJ 



xjg(Qhx,o) 



(13) 



This equation can be derived by considering a linear re- 
sponse on a time-dependent magnetic field that is adi- 
abatically switched on at t = —00, and then suddenly 
switched off at t = 0, i. e. B(i) = e St 6(-t)B with 5 -> 



S x (p,t) = S± (p,0)V(Q hx ,t), 



(10) 



(see, e. g., Ref. [33j for similar calculations). Usually the 
SO coupling is weak on the Fermi energy scale, which 
implies Qhx kp, where kp is the Fermi momentum. 
Therefore in most situations one can safely replace the 
static uj = response function in Eq. (|f 3[) by the macro- 

r en 

scopic Pauli spin susceptibility, Xpp{Qii*,fy = xp- 

Equations (fT"2|) and ([T3| are the result of the spin- 
density dynamics duality and give the exact evolution 
of the uniform spin density. The problem is solved 
by mapping the spin relaxation in the physical system 
to the "washing out" an inhomogeneous spin texture 
in a dual system without SO coupling. The real spin 
relaxes because of SO-induced precession and random- 
ness introduced by disorder, phonons, and interelectron 
interactions. [HI For the transformed spin, it is the evolu- 
tion of the nonuniform spin density distributions, again in 
the presence of the disorder and interactions between the 
carriers. An interesting exact feature of the pure gauge 
SO coupling (in addition to the well known anisotropy) 
is the absence of the spin precession - the vector Sj_(f) 
is always collinear to its initial direction. In the trans- 
formed picture this is related to the diagonal structure of 
the spin response in a nonmagnetic electron gas. For the 
real system this translates to the fact that spins of elec- 
trons with opposite momenta precess around the h-axis 
in the opposite directions with the same rate. 
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At the level of the random phase approximation the 

rol 

spin response function Xpp{Q> w ) is equal to the den- 
sity response function x(Qi w ) of a noninteracting, but 
possibly disordered and/or coupled to phonons electron 
gas, while the Pauli susceptibility xp is proportional 
to the compressibility dn/dfj., with n and /z being the 
electron concentration and chemical potential, respec- 
tively. Hence the spin diffusion Green's function entering 
Eq. (TT2"]) reduces to 



T>(Q hx ,uj) = — 

1UJ 



x(Qhx,^) 

dn/dfj, 



1 



(14) 



which is exactly the density diffusion Green's function. 
Therefore in this physically important case the spin re- 
laxation is mapped to the ordinary density diffusion. 

Now we apply Eqs. ([T4"|) and (TT2"]) to a noninteracting 
disordered 2D electron gas with a momentum relaxation 
time t, and study possible regimes of spin dynamics. The 
density-density correlator x(Qhx,u) can be obtained ei- 
ther diagrammatically or by solving the kinetic equation. 
In the semiclassical regime, corresponding to the summa- 
tion of ladder diagrams, one obtains: 



V(Q hx ,uj) = 



1 - K. (Qhx, to) ' 



£(Qhx,w) = /- r 

2tt J l-i 



dO 



\LjJT + ifL n T cos ( 



(15) 
(16) 



where the only SO-dependent parameter in the problem 
£l so = QhxVF {vf is the Fermi velocity) is the maximum 
spin precession rate, and Q S oT = £Qh x (electron mean 
free path £ = vft) characterizes the relaxation regime. 

We begin with f2 so T <C 1 regime, which coresponds to 
a pure diffusion, studied in the coordinate representation 
in Ref . [ljj ■ Here the diffusion Green's function, Eq. (fTo]) , 
reduces to 



(Qhx,w) 



DQ 



2 

hx 



(17) 



where D ~ Vpr/2 is the diffusion coefficient. Inserting 
D (<2hx, oj) of Eq. (TTT]) into Eq. (TT5J we obtain: 

(18) 



S ± (t) = S ± (Q)cxp(-DQlj), 



which exactly corresponds to the Dyakonov-Perel' mech- 
anism with the spin relaxation rate T s = DQ^ X . More- 
over, the factor 1/2 in the definition of D acquires an 
interesting physical meaning in terms of the spin pre- 
cession: it corresponds to the angular averaging of the 
precession rate (0^ (k)) = f2j? /2. 

The opposite, clean limit Q S oT 3> 1, in terms of the 
dual (transformed) system corresponds to a reversible, 
purely ballistic washing out the helix texture. In this 
regime V(Qh^,uj) w /C(Qhx,^) (see Ref.[35|) and the 
integration in Eqs. (fT2|) and (flB")) yields: 

S± (f) = Sx(0) Jo (n so t) = Sx(0) Jo (Qhxvpt) , (19) 
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FIG. 1: (Color online.) Time dependence of the spin for 
different parameters of SO coupling, shown near the plots, 
with s(t) defined as Sx (t) = s(t)8± (0) /5x(0). 



where Jo (f2 S oi) is the Bessel function. The same result 
can be derived directly from the microscopic spin preces- 
sion with a k-dependent rate fl so (k) = f2 so cos <f> 1 where 
4> is the angle between k and v. Indeed, the net result of 
the inhomogencous precession reproduces Eq. (|19p . 



Sx it) = Sx(0) / cos(fi so icos 



2tt 



Sx(0)J (fi so i). 

(20) 

In contrast, in the systems with pure Rashba or Dressel- 
haus coupling, spin precession rate does not depend on 
the direction of momentum. As a result, in the ballistic 
regime, the z-component of the total spin demonstrates 
cosine-like rather than the Bessel function time depen- 
dence. 

An intermediate regime of fl so r ~ 1, can be investi- 
gated numerically. The results are presented in Fig[T] for 
different parameters. One can clearly see a crossover from 
the oscillating Bessel function-like behavior to the expo- 
nential Dyakonov-Perel' decay. At short times, the be- 
havior of spin is universal: Sx (t) = Sx (0) (l — fig Q i 2 /2) 
due to the unperturbed precession of the spins. For the 
density dynamics the universal short-time behavior is a 
direct consequence of the /-sum rule,[36| as can be seen 
by expanding Eq. (Tl"2"|) at t — > 0. 

Analysis of Eqs. ([T2" ]) .(fl~5 ]) . and (fTBj) shows that f2 so r = 
1 is a critical point. With the decrease in f2 so T in the 
range il so T > 1, the first zero of Sx (t) rapidly shifts to 
larger times, with negative value regions becoming very 
shallow. At ft so T < 1 zeroes of Sx (t) disappear and the 
dynamics is a pure decay. 

The gauge transformation approach allows to analyze 
systems, where a direct treatment of the SO coupling 
would cause difficulties. The first effect we consider is 
the influence of the orbital motion in a nonquantizing 
magnetic field B along the z-axis on the spin dynamics. 
We assume that due to a small g— factor of electron, the 
Zeeman coupling to the magnetic field does not cause a 
relevant spin precession. If fl so r <C 1, the density evolu- 
tion at B = is diffusive, and the electron mobility de- 
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creases with B due to the Lorentz force as (l + w 2 t 2 ) , 
where uj c = \e\B/mc is the cyclotron frequency. By 
the Einstein relation, the diffusion coefficient is renor- 
malized by the same factor D(B) = D(0)/ (l + lj' 2 t 2 ). 
Hence the spin relaxation rate in Eq. (fT8|) decreases as 
T S (B) = r s (0)/ (l + w 2 t 2 ), which reproduces the results 
of the direct quantum kinetic theory. [33 ] For illustration, 
we consider the limit uj c t 1 at short times t <C t in 
a more detail. Here the trajectories of electrons are very 
close to circles and the spin-independent kernel in Eq. (|9"j) 
can be represented as (cf. Ref. 351 ]) 



V(p'-p,t) 



1 



2nd(t) 



S(\p'-p\-d(t)), (21) 



with the displacement d(t) = 2R C |sin (o; c i/2)| , where 
R c = vp/uj c is the cyclotron radius, and R c Qhx = 
il so /cj c . Straightforward integration in Eq. (|9|) yields (cf. 
Eq.®) 



S± (*) = S±(0) Jo (Qhxd(t)) : 



(22) 



where Qhxd(t) = 2Q so |sin(cj c t/2)| /w c . For a very weak 
field (in a very clean system) w c <C £l S o, the result re- 
produces Eq. (fT9|) . as expected. In the opposite limit 

In terms of spin pre- 
as very fast 



3> f! so , no relaxation occurs. 
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cession this can be understood [2 z 
changes at the frequency uj c in the direction of the SO 
field, keeping the total spin out of relaxation. In terms 
of the nonuniform density dynamics, this implies that 
the electrons, forced to circulate around small radius cy- 
clotron orbits, can not spread out to destroy large scale 
~ 1/Qhx 3> R c density variations. At long times, the dif- 
fusion behavior takes over, and the relaxation becomes 
exponential. 

As a second example we briefly discuss the effect of 
weak localization on the spin relaxation 4^, 41 1 by con- 
sidering w-dependent renormalization of the momentum 
relaxation rate with the correction: (42 - 44 1 



M = - 



i 



d 2 q 



1 



Tdn/d^J (2ir) 2 Dq 2 — iw 



(23) 



This correction arises due to the enhanced return prob- 
ability, which slows down of the density dynamics and 
eventually leads to the algebraic tail in spin relaxation at 
long time as S_l (t) ~ 1/t. [Il[ In terms of spin preces- 
sion, the enhanced backscattering slows down the spin 
relaxation because upon the return to the initial point 
the direction of the electron spin remains the same. 

To conclude, we have shown that in a wide class of 
systems the non-Abelian gauge field description of SO 
coupling reveals the duality of experimental observables 
and ensures the exact mapping of the spin dynamics 
to the density evolution. The evolution is described in 
terms of the responce to an external perturbation with 
the wavevector equal to the spin helix wavevector Qi lx . 



We presented explicit results for a weak SO coupling with 
Qhx <C kp valid for all systems of interest (this restriction 
is not required in general). This exact mapping opens a 
venue for understanding the whole class of practically im- 
portant systems through better studied, and more simple 
properties of the systems without SO coupling. 
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